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Abstract 

We calculate Aa, T,a and Ha potentials from the nuclear-matter G-matrices of 
the SUq quark- model baryon-baryon interaction. The a-cluster wave function is as- 
sumed to be a simple harmonic-oscillator shell-model wave function. A new method 
is proposed to derive the direct and knock-on terms of the interaction Born ker- 
nel from the hyperon-nucleon G-matrices, with explicit treatments of the nonlo- 
cality and the center-of-mass motion between the hyperon and a. We find that 
the SUq quark-model baryon-baryon interactions, FSS and fss2, yield a reasonable 
bound-state energy for ^He, —3.18 ~ —3.62 MeV, in spite of the fact that they 
give relatively large depths for the A single-particle potentials, 46 ~ 48 MeV, in 
symmetric nuclear matter. An equivalent local potential derived from the Wigner 
transform of the nonlocal Aa kernel shows a strong energy dependence for the in- 
cident A-particle, indicating the importance of the strangeness-exchange process in 
the original hyperon-nucleon interaction. For the Sa and Ha potentials, we only dis- 
cuss the zero-momentum Wigner transform of the interaction kernels, since these 
interactions turn out to be repulsive when the two isospin contributions for the 
TiN and interactions are added up. These components show a strong isospin 
dependence: They are attractive in the isospin 1 = 1/2 (Sa) and 1 = (Ha) com- 
ponents and repulsive in / = 3/2 (Sa) and 1 = 1 (Ha) components, which indicate 
that £ and H potentials could be attractive in some particular systems such as the 
well-known gHe system. 
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1 Introduction 



Interactions between the octet baryons (B 8 = N, A, E and S) and the a 
cluster are important ingredients to consider the possible existence for various 
kinds of light hypernuclei through detailed microscopic cluster-model calcu- 
lations. If reliable effective B$N interactions were known, one could easily 
calculate the B 8 a potentials, using the standard cluster-model techniques. 
Unfortunately, this is not the case except for B 8 = N, since the bare B 8 N in- 
teraction itself is not well known especially for the E and S hyperons, because 
of the technical difficulties of strangeness experiments. Even if these bare in- 
teractions were eventually known, we further need to develop the procedure 
to link the bare interactions and effective interactions through some effective 
interaction theory such as the G-matrix formalism. 

We have recently developed the QCD-inspired spin-flavor SUq quark model 
for the baryon-baryon interaction [1], which is a unified model for the full 
octet-baryons [2], and have achieved accurate descriptions of the NN and YN 
interactions [3]. In particular, the NN interaction of the most recent model 
fss2 is accurate enough to compare with modern realistic meson-exchange 
models. These quark-model interactions were used for the detailed study of 
few-baryon systems such as 3 H [4,5] and \H [6], and also of some typical A- 
hypernuclei, \Be [7,8] and A ^He [9], through a newly developed three-cluster 
Faddeev formalism [10,11] and G-matrix calculations [12,13,14]. We can now 
use these baryon-baryon interactions to calculate not only the Aa interaction, 
but also Ea and Ha interactions, assuming the harmonic-oscillator (h.o.) shell- 
model wave function for the a-cluster. 

There are, in fact, a couple of advanced procedure to derive the interactions 
between a single baryon and finite nuclei based on the density-dependent 
Hartree-Fock theory of G-matrix interactions [15,16,17]. These approaches, 
however, use the localized G-matrix interaction in the configuration space 
and the center-of-mass (cm.) coordinate system connected to the target nu- 
cleus. For the B 8 a interactions, the cm. correction is quite important. In this 
paper, we will derive B 8 a Born kernels, directly starting from the G-matrix 
calculation of the quark-model baryon-baryon interactions. We deal with the 
nonlocality of the G-matrix and the cm. motion exactly, using the cluster- 
model approach, at the expense of the self-consistency between the a-cluster 
formation and the G-matrix interaction. The G-matrix is pre-determined by 
solving the Bethe-Goldstone equation in symmetric nuclear matter, and the 
momentum-dependent single-particle (s.p.) potentials as well [12]. The Fermi- 
momentum is assumed to be the standard value k F = 1.35 fin -1 , but the 
obtained B 8 a interactions depend on this choice rather monotonously, except 
for some special cases such as the Aa LS interaction. The treatment of the 
starting-energy dependence and the non-Galilean invariant momentum depen- 
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dence of the G-matrix is explicitly discussed. We believe that these approxi- 
mations are good enough to understand the unknown £a and Ha interactions 
qualitatively, starting from the quark-model predictions of various baryon- 
baryon interactions. For the Aa interaction, we compare the predictions by the 
present approach with some available phenomeno logical Aa potentials [7,18], 
obtained by various methods. Another application of the present approach to 
the Na interaction will be published in a separate paper, since this system 
involves an extra nucleon-exchange term. 

Since the obtained B 8 a interactions are all nonlocal, we calculate the Wigner 
transform from the B 8 a Born kernels. We find that the momentum dependence 
of the Wigner transform is very strong, and the procedure to find effective lo- 
cal potential by solving the transcendental equation is necessary to obtain a 
local-potential image for the Aa interaction. The £a and Ha interactions are 
repulsive, although the isospin 1 = 1/2 (for Ha) and 1 = (for Ha) contribu- 
tions of the HN and EN interactions are attractive. As the first step to study 
realistic Ha and Ha interactions, we will only discuss the zero-momentum 
Wigner transform in this paper. Although these B 8 a interactions are complex 
due to the imaginary part of the underlying G-matrices, we discuss only the 
real part. The spin-orbit B 8 a potentials are naturally obtained from the LS 
and LS^ components of the G-matrix invariant interaction. 

The central Aa, Ha and Ha potentials were calculated by many authors 
with other (usually more crude or purely phenomenological) approaches. The 
Aa potential was calculated numerously, and different shapes were consid- 
ered. For instance, see Refs. [19,20,21,22]. Some phenomenological Ha and 
Ha potentials are found in Ref. [23,24,25,26]. The comparison of our results 
with these potentials are, however, not easy, since our Wigner transform is 
momentum-dependent and the solutions of the transcendental equations are 
strongly energy- dependent. We will only point out if some resemblance be- 
tween them are found. The importance of the isospin dependence of the HN 
and HA" interactions to the S-nucleus and H-nucleus potentials has been dis- 
cussed by many authors from the phenomenological aspects, for example, in 
Refs. [27,23,25], and by Rijken and Yamamoto from the viewpoint of G-matrix 
calculations in Ref. [28]. 

The organization of this paper is as follows. In the next formulation section, 
we first give in Section 2.1 the basic folding formula for the a-cluster based 
on the separation of the B 8 a Born kernel to the spin-isospin factors and the 
spatial part. The treatment of G-matrix variables, the starting energy and 
the cm. momentum, will be carefully discussed. A convenient transformation 
formula for the rearrangement of relative momenta is given in Section 2.2, by 
which the partial wave component of the B 8 a Born kernel is explicitly given 
both for the central component and the LS component. The folding formula in 
the partial-wave expansion and the partial-wave components of the B 8 a Born 
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kernel are explicitly given in Section 2.3. A procedure to calculate the Wigner 
transform is given in Section 2.4 with a couple of approximations. One of the 
approximations for the LS component yields a simple factor for the strength of 
the LS potential, which corresponds to the well-known Scheerbaum factor S B 
[29] in nuclear matter. The Section 3 is devoted to the results and discussion; 
first in Section 3.1 the Aa central and LS potentials both in the T-matrix 
approach and in the Wigner transform approach. The isospin dependence of 
the Set and Ea potentials is discussed in Section 3.2. Section 4 is devoted to 
a summary. The invariant G-matrix for the most general B 8 B 8 interaction is 
discussed in Appendix A. 



2 Formulation 



2.1 a -cluster folding for the G -matrix invariant interaction 



The B 8 a Born kernel for the (Os) h.o. a-cluster wave function is calculated 
from 

V(q„ Qi ) = ( 5{X G )e^f r XB K | £ | 1 • e^ r X bK > , (2.1) 

i=2 

where xb is the spin-isospin wave function of B, <p a the internal wave function 
of a, and G\j the BN G-matrix acting on the particle % — 1 (B) and j '• = 2 
- 5 (nucleons). We use the short-hand notation B to specify one of the octet 
baryons, B 8 = N, A, S or S. In Eq. (2.1), it is important to calculate the 
Born kernel in the total cm. system by inserting 5(Xc) and 1 in the bra and 
ket sides, respectively [30] , since the G- matrix interaction G\j is non-Galilean 
invariant. Namely, the two-particle G-matrix which satisfies the translational 
invariance is parametrized by 

(Pi,P2\G\p[,p' 2 ) = 5(K - K')-^G(p,p';K,u;,k F ) . (2.2) 

Here, ui is the starting energy, K = \K\ is the magnitude of the cm. momen- 
tum, kp is the Fermi momentum of nuclear matter, and the relative momentum 
p (and also p' etc. with primes), is given by 



n £Pl ~ P2 

K = Pl +p 2 , 



Pi 



P2 



K + p , 



K-p , 



(2.3) 



4 



with £ = {Mn /Mb) < 1 being the mass ratio between the nucleon and the 
baryon B\^j We assume a constant kp = 1-35 fm -1 in this paper unless oth- 
erwise specified, so that we will omit this index in the following. In fact, 
the G-matrix in Eq. (2.1) contains the exchange term. It is, therefore, conve- 
nient to use the isospin sum of the invariant G- matrix, G BB (p, p'\ K,lu) with 
I — Ib + 1/2 and Is — 1/2 (Ib is the isospin of B), defined through 



G bb (p,p';K,uj) 

= ( [BN\ IIz | G(p, p'; K, u) - G(p, -p', K, u) P a P F | [BN\ IIz ) 



'1 + 5 



B,N 



9o + aLi* 7 ! ■ cr 2) +h^in- (<ri + cr 2 ) + hi in 



<?2 



+ 



(2.4) 



Here n = [p'xp]/ (p'p sin 9), and the invariant functions g$ (central), g^ s (spin- 
spin), h* Q (LS), hi (LS^), etc. are functions of p 2 , p' 2 , cos 6 = (p ■ p'), K, 
u, and k F . These are expressed by the partial-wave components of the BN 
G-matrix as (see Appendix D of Ref. [31]) 



1 



J£S 



>Gs J i,st p i( cos( 



1 ^ (2J+1) 

A 2-~i 



j J(J + 1) 



h 1 



i[25(5 + l)-3] ^ 
Gilu P}(cos9) 
u+lu+i P}+i(cosO) - (J+^GH^j^ PU(cos( 

P}(cos9) , 



+J G IJ 



1 (2J+1 

A 2-^1 



J(J+1) 



i j 

0J.1J 



(2.5) 



where the argument p, p', K, u, kp and subscripts for the baryon chan- 
nels are omitted for the typographical reason. In the LS term, P](cos9) = 
(sin 9) P'j (cos 9) with J > 1 is the associated Legendre function of the first 
rank. The invariant G-matrix for the most general B$Bs interaction is dis- 
cussed in Appendix A. 

In order to calculate the spin-isospin factors, it is convenient to introduce an 
isospin pseudo-exchange operator of the BN system by 



Pr 



2/b + 1 



[1 + t b - r N ) 



(2.6) 



1 Note that ( used in Appendix B of Ref. [7] is the inverse of £. 
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with the isospin matrix elements 

2I B 

(tb ■ t n ) = < 

k -2(1 B + 1) 
and write the invariant G-matrix as 

G(p, p'; K, u) = G J B = /B+1/2 (p, p'; K, u ) 
We need to calculate 



for 



/ = < 



Ib + 1/2 
Ib ~ 1/2 



(2.7) 



— — +G BB B (p,p;K,uj) — - — 



j=2 V 



1 ±P T 



\XBXa) 



(2.8) 



(2.9) 



with the spin factors u^f 



central 



LS(-) 
12 



i, ^ 



(<ri • cr 2 ), c^ 5 = cri + <x 2 , and 



(Ti — <x 2 . In Eq. (2.9), Xa is the spin-isospin wave function of a; i.e., 



0a = Xa0a pace - Then we find that non-zero matrix elements are 
21 + 1 



/ central 



= 2- 



2/ f 



y-/L5 _ y/LS<-) _ q ^+ ] 

B B 2J B + 1 



CTi . 



(2.10) 



On the other hand, the spatial part is calculated in Appendix (B6) of Ref. [7]. 
It is convenient to express this formula as 



V space (q f , q { ) = { 5(X G )e tq f r <f) s v ace | G spacc | 1 • e l ^- r 0^ pace ) 

xG^(p+±k,p-±k;(l + t)\q-p\,u>} , (2.11) 

where 

k = q f - qi , q=^(q f + q l ) , (2.12) 

are the momentum transfer and the local momentum of the B 8 a system. By 
assuming /i^in and hi. in for G spacc , we finally obtain 
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gi(p+\k,p-\k-(l + ()\q-p\,u) 1 ^ i3 ^ 

^ (p + \k, p - ife; (1 + £)\q - p|, u) i[p x fc] • o-i ^ 

where h 1 = Hq + h T _ and [p x fc] = [p x fc]/| [p x fc] |. 

It should be noted that the cm. momentum of the two interacting particles, 
K = (l+£)(qr — p), in Eqs. (2.11) and (2.13) implies that the local momentum 
q = {Qf + <?i)/2 in Eq. (2.12) plays the role of the incident momentum of the 
first baryon B from Eq. (2.3), and p the local momentum of the two-particle 
system, which is now an integral variable. This is a consequence of fixing the 
cm. motion of the B 8 a system as in Eq. (2.1), and is a special situation of 
the direct and knock-on (exchange) terms. The G-matrix depends only on 
the magnitude K = \K\, since we make an angular average in the G-matrix 
calculation [12,13]. The G-matrix value is, therefore, specified by K and u, or 
alternatively by the incident momentum qi = |q 1 | and the relative momentum 
q = \q\ between B and iV; i.e., G(p ) p']K ) uo) = G(q + k/2,q — k/2;q 1 ,q) with 
k = p — p' and q = (p + p')/2. If q\ and q are specified, K is determined by 
the angular averaging, and q2 is determined from q 1 and K. Then the starting 
energy uo is determined as the sum of the relative kinetic energy of two particles 
and the s.p. potentials for B 8 and iV at q\ and q 2 , respectively. We therefore 
choose qi = \q$ + qfJ/2 and q = \p\ in the G-matrix in Eqs. (2.11) and (2.13). 
In order to carry out this rather involved calculation, we will develop in the 
next subsection some kind of transformation formula for the rearrangement of 
relative momenta in the partial-wave components of nonlocal kernels. 



2.2 A transformation formula of the nonlocal kernel in the momentum rep- 
resentation 



The folding formula in Eq. (2.11) implies that expressing the G-matrix inter- 
action G c (p, p') and G LS (p,p') i[p' x p] • S with the subsidiary momentum 
variables 



k=p-p' , 

q = ^(p + p) , 



p = q + -k , 



p 



q~- 2 k , 



(2.14) 
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is convenient for the a-cluster folding. We express these kernels using the 
calligraphic letters; i.e., 



g C \k,q)=G C (p,p>) , 

Q LS (k, q) t[q~7k] ■ S = G LS (p, p') i[p^Tp] ■ s 



(2.15) 



By using this notation and the qi, q notation for K, u, discussed in the last 
subsection, G s P ace in Eq. (2.11), for example, becomes Q spacc (k,p;qi,p) with 
Qi = \Qf + 9il/2- In the following, we omit the argument K, uj or q 1 , q for 
simplicity, unless the dependence becomes crucial for the relationship under 
consideration. After making partial wave decomposition in Eq. (2.15), we can 
easily carry out the integral over p. The resultant expression is denoted by 
V c {k, q) and by V LS (k, q)i[q x k] ■ S. The desired Born kernels, V c (q f ,q i ) 

and V LS (qj, q^ i\q i x qf]-S, are obtained from another transformation related 
to Eq. (2.12). 

Our task is, therefore, to relate the partial-wave components G\(k,q) and 
Gf(p,p'), and also G\ S (k,q) and Gf s (p,p'), which are defined through 



g c (k,q) = Y,(^ + l)Qx(k,q)P x (k-q) , 

A=0 

oo 

G c \p,p') = £(2€ + l)Gf(p,p')P e (p ■ p') , 
g LS (k, q) = £) (2A + l)G L x s {k, g)P A 1 (fc ■ q) , 

A=l 

oo 

G LS (p,p>) = ^+1)G^ S (p,p')PKp-p') 
e=i 

We generalize the transformation Eq. (2.14) as 




/ 



a (3 \ I k 



with a5 — /?7 = 1. A simple calculation gives 

Ox(k, q) = flW + Vll C dx 9jMl Px ( x ) ge (k, q- x) , 

where ge(k, q; x) = (pp') e Pe(p • p') with 

p = \J a 2 k 2 + f3 2 q 2 + 2af3kqx , p' = \J^ 2 k 2 + 5 2 q 2 + 2^5kqx 



(2.16) 



(2.17) 



(2.18) 



(2.19) 
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The transformation of the LS components is similarly carried out by using 
the orthogonality relations of P}{x) and 



p/(cos0)= E (2^ + l)/V(cos0) , 

VS1Iiy/ <'=<-l/-3,-,l or 

Px-i(x) - p x+1 (x) = X(XTi) v/r3 ^ p fa) ■ ( 2 - 2 °) 

For the numerical integration over x in Eq. (2.18) etc., it is convenient to use 
the spline interpolation 

<%(p,lO = E s i(p) s iV)G?(Pi,Pi) » ( 2 - 21 ) 

since the G-matrix calculation itself is very much time consuming. After all, 
we obtain the following formula for the transformation of the partial-wave 
components: 



0?(M) = E E (2€ + i)i$*<3W;) (q = c,ls) 

i,j £=0 or 1 

p«« - _J^_ V (2/ +1)1/" <ix M 



r=^-l,^-3,- ,1 or 

x [Px-i(x) - P x+1 (x)\ g e (k, q; x) , (2.22) 



where 



g £ (k,q;x) = (pp') 1 l P t (p ■ p') 

= y y t-iY'+ti' ( 2£+ x ) ! a /i^ 2 yi'^' 

t^=th'ttj=i ! (2£ 2 ) ! (27?) ! (24') ' 



x tf 1+tl '<^ +t '' E(40£i'0|£'0) (404'0|/0) 



CI C2 t 

4' 4' ^ 



iM^) • (2-23) 



The summation over £ in Eq. (2.22) is from I = for Q = C and from i = 1 
for n = LS. 

It is important to note some symmetries possessed by G\(k, q). From the time- 
reversal symmetry, the G-matrix G n (p,p') is symmetric for the interchange 
of p and p'. Since this interchange corresponds to k — > — fc, the transformed 
partial-wave components, G\(k,q) and G\ S (k,q), are no-zero only for A = 
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, 2 , 4 , • • • and A = l,3,5,---, respectively. For the coefficients in Eq. (2.17), 
we can show gi(2q,k/2;x) = (—l)g^(k,q;x) in Eq. (2.23). This property im- 
plies that, if G^(p,p') is transformed to Qjf (k, q), then (— l) e Gf(p,p') is trans- 
formed to Q^(2q,k/2). In the later application of the present formalism to 
na resonating-group method (RGM), we will find that the knock-on term is 
obtained by simply replacing G\(k,q) to Q^(2q,k/2) in the direct term. In 
fact, the knock-on term is already included even in the hyperon a system as 
in Eq. (2.4), which is the strangeness exchange term of the hyperon-nucleon 
interaction. In other words, the direct and knock-on terms are treated on the 
equal footing in the present formalism to deal with the invariant G-matrix 
interaction. For the transformation from V n (k, q) to V n (qj, q^, we find that 
the latter Born kernel is symmetric with respect to the interchange of q$ and 
qi, which is the consequence of Eqs. (2.19) and (2.22) for the transformation 
coefficients (see Eq. (2.12)) 



(2.24) 



2.3 Partial-wave expansion of the B 8 a Born kernel 

If we use partial- wave decomposition in Eq. (2.16) and the similar expansions 
for V n (k, q) = V n (qf, qj, the folding formula Eq. (2.11) becomes very simple. 
For both of the VL = C and LS terms, it is given by 




v?(M)^-^-|(i + ,)V}^ + « 

r 2, f 2(1 +o 2 2 i . /(i +40 

x ^ p 2 dp exp I ^ p 2 \ i x I ^ ^- ^ 



where i\(x) = i x j\(—ix) is the spherical Bessel function of the imaginary 
argument. For the proof of the LS'-term folding, we use a simple formula 




P}(k-q)i[^] v = (-l)'^^^. [Y e (k)Y e (q)] iu , (2.26) 

which is derived by using 

P} +1 {x) - P}_ x (x) = (21 + 1) VT^P e (x) , (2.27) 

and 
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P t (k -q)i[qx k] u = (-1 



e Air kq 



1 (£ + !)(£ + 2) 



x 



Y e+1 {q)Y e+1 (k) 



V3(2£+ 1) IV 2/ + 3 



'(I -1)1 
2£-l 



(2.28) 



The final step to derive the partial-wave components Vp{qf, qi) in 
V(q f , ft) = V C (q f , ft) + V LS (q f , ft) ifoTj/] ■ S , 

oo 

V C (q f , ft) = E( 2 ^ + 1) Vf(q,, q % ) P e (q f ■ ft) , 



oo 



V LS (q f , ft) = £(2£ + 1) V £ LS (q f , q t ) P}(q f ■ ft) , 



e=i 



(2.29) 



is carried out by using the transformation formula in the preceding subsection 
with the coefficients in Eq. (2.24). In the jj-coupling scheme, we also need 
V/(qf,qi) in the partial-wave expansion 

V(q f , ft) = 4tt £ V/(q f , ft ) ^ y { q )JM (<Z/5 SP^) ^i )JM (ft; spin) , 
je m 

(2.30) 

which is given by 



V e J {q f , qi) = V e C (q f , q t ) + V t Lli (q f , Qi ) (L ■ S) je , 



LS, 



{L-S) Jt = 



J(J + !)-£(£+!) 



2 

e±i 

2 



for J 



+ 1 
_ i 

2 



(2.31) 



In Eq. (2.30), y^ JM (ft - spin) = [Yg(ft) Xi]jM is the angular-spin wave func- 
tion of the B s a system. For the proof, we again use Eq. (2.26). 

In summary, the Born kernel of the B 8 a system is obtained by using Eqs. 
(2.22) and (2.25), starting from 
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= \ g (I£t) (IS) = o. i. V • •) , 



(C=l, 2, 3, •■) . (2.32) 



For B s = N, we should multiply the factor 2 to include the knock-on term. 
We first use Eq. (2.22) with the coefficients Eq. (2.17) and transform the above 
Gp(p,f/) to G\(k, q). Then, the a-cluster folding by Eq. (2.25) yields Vj?(fc, q) 
from G\{k, q). The second transformation from V x (k, q) to V x (qf, qi) is carried 
out with the coefficient Eq. (2.24). The selection of q\ and q in G n (p,p'; q±, q) is 
now almost apparent. We choose q± as q in the folding formula Eq. (2.25). The 
relative momentum q is actually p in Eq. (2.25). Therefore, the nest structure 
V$(h,q) D G%(k,p; qi = q,q = p) and G%(k,q;q 1: q) D Gf(p,p';q u q) is incor- 
porated in the computer code. Namely, we first specify k and q in Eq. (2.25). 
Then G\(k,p;qi = q,q = p) is generated from G\(k,q;qi,q) for a general q, 
which is obtained from the transformation formula, Eq. (2.22), by using the 
complete off-shell G-matrix Gf(pi,pj; qi, q). 

2.4 Wigner transform 

The present formalism is convenient to calculate the Wigner transform V^(r, q), 
since they are essentially the Fourier transform of V n (k,q). We define these 
through 



Note that the LS term is defined by L w • S = [rxq]-S, instead of i[r x q]- S. 
If we apply the partial-wave expansion 




= j^fl dk etk ' r { vC(fcj q) + yL5(fc ' q) i[q ^ ] ■ s ) 

= V^(r,q) + V^ s (r,q)[rxq].S . 



(2.33) 



oc 



V^(r,q) = j:(2\+l)V^ 



(r, q)Px(cosip) 



A=0 



V&>,g) = £(2A+l)V&i(r, g ) 



A=l 



oc 



( 



sin ip 



1 



) 



P x (cos ip) 



(2.34) 
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with cosy? = (f • q), we can easily derive the partial- wave components as 



A-TT roo 

V wx(r,q) = j^y 3 J Q k 2 dk * A jx(kr)V°(k, q) , 

V ^ q) = VrJ^fL k2dkiX ~ l h{kr)V L x S {Kq) . (2.35) 

Here we use Eq. (2.26) again for the derivation of the LS term. Since V^ x (r, q) ^ 
for A = 0, 2, 4, • • • and V^f (r, q) ^ for A = 1, 3, 5, • • • (see Subsec. 2.2), 
the A = and A = 1 terms become the leading terms for the central and LS 
Wigner transform, respectively. In fact, in the q = case, we find that only 
these leading terms survive for the zero-momentum Wigner transform. It is, 
therefore, a good approximation to retain only A = and A = 1 terms in 
Eq. (2.34): 



Ajr roo 

VC(r,q)~-— k 2 dkj (kr)V^k,q) , 
(Ztt)- 3 jo 

V w(r,q)~j^j- 3 J o k 2 dk- 3l (kr)-^ s (k,q) , (2.36) 
which we use throughout this paper. 

The zero-momentum Wigner transform is a convenient tool to estimate the 
strength of the interaction. From the folding formula Eq. (2.25), we can cal- 
culate V$(k, 0) and (3/q)Vi S (k, q)\ q =Q. This process yields 



V#(r,0) = - [°°k 2 dk e-^ k2 3o (kr) 

7T JO 

/2(i + o 2 \^ r° 2j m+ov - c ., , 

Jo 



37TZ/ 



3f 7r Jo r 



(2.37) 



For the LS component, it is easy to introduce another approximation, which 
gives a simple factor for the LS strength, similar to the Scheerbaum's factor 
[29]. For this purpose, we express Gi S (k,q) in the original form 
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(*, q) = % f(2i +1)1 f d X [1 - P 2 (x)] 

o z i pp 

Pi {cose) , 



(-L- 

Vsin# 



(2.38) 



where cos# = (p ■ p) with p = q + fe/2, p' = q — k/2 and x — (k ■ q). We 
neglect the k dependence except for the front factor kq/3 and set k = 0. Then, 
p, p' — > g, cos# — > 1, and the x integral can be performed. By further using 
p;(l) = £(£ + i)/2, we obtain 



6 i=i q z 



(2.39) 



If we set 



1 



<?(?) = 4 £€(€+ 1)(2€+1) g) 



£=1 

1 ^ / 27 + 1 



E 



2/r + 



T j (2 J + 1) {-Gft,! q) — (J + 2)G[ J J+1<1 J+1 (q, q) 



l)G^_i,u_i(?,9) + //(J+l) 
27+ 1 



2E 

u 



2/b + 1 



K(p,p') + hUp,p') 



G[ J Jfi Aq,q) + G^ 1J (q,q) 
1 



sin# 



(2.40) 



we can write §i S {k, q) ~ (2k/3q)G(q). Using this and the integration formula 
2 r 00 , , „ 1 



- / k 3 dk e sL fc2 -jAkr) = —= ( — ) ' e 

7T JO r y' 7T V 3 / 



(2.41) 



we find 



^o,^(l±i^±iilA (f) %^ 



X 



'2(1 + 0' 



2 i 2(1+0 a 2 ^/ \ 

g 2 <ig e « G(g) 



y 3ttv J jo 

We can write Eq. (2.42) in the way similar to the Scheerbaum's formula: 



(2.42) 



TT( \ 71 Q 1 ^( r ) 

2 r dr 



(2.43) 
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For the B 8 a system, the density p(r) of the a-cluster is calculated as 

P(r) = (0 a |E5(r-^)|0 Q )=4g)%-|- 2 , (2.44) 

where ^ i = Xi — X a is the coordinate of the nucleon i, measured from the 
cm. of the a-cluster. Then, we find that the integral Eq. (2.41) is nothing but 
(—n)(l/r)(dp(r)/dr). Thus Eq. (2.43) becomes 

U(r) = S B *= (|) S . (2.45) 

Similarly, we can write Eq. (2.42) as 

V™(r, )=S B -^(^f e^' , (2 .46) 

and obtain 

This expression corresponds to Eq. (50) of Ref. [13], which is the Scheerbaum 
factor of finite nuclei derived from G-matrix calculations: 

1 £ 3 o rSmax 

Sb <*> = 2^171 (M5< 1+f) I d * W ^M . (2.48) 
The difference between Eqs. (2.47) and (2.48) is the weight function, 

W(q) = q 2 e~ 2J1 ^ Lq2 < ► W(q u q) , (2.49) 

and an extra front factor in Eq. (2.47). This enhancement factor appears, 
since in the standard G-matrix calculation of single-particle potentials the 
cm. motion of the total system is not correctly treated. In the B^ot system, 
this approximation affects the result appreciably, which is discussed in Ref. [8] . 
The two weight functions in Eq. (2.49) may seem to be fairly different, since 
W(q 1: q) with q 1 = is given by W(0, q) = 9(q — k F /(l +0)- However, this is 
not the case, since in the small q region (the low-energy region) G(q) is anyway 
very small. We can further approximate Eqs. (2.47) and (2.48), by calculating 

Q = \J (Q 2 ) with each weight function, and by replacing q in G(q) with this q. 
The average momentum q is given by 

' = ¥(lTf) for w(q) ' « = ^(lT?) for ,1/(0 ' 9) ' 

(2.50) 
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For W(q) with v = 0.257 fin -1 , q ~ 0.38 fm -1 , and for W(0, with k F = 
1.2 fm -1 , q ~ 0.34 fm -1 . These values are very close to q = k/2 = 0.35 fm -1 , 
calculated from the Scheerbaum's estimation for an average momentum trans- 
fer k ~ 0.7 fm -1 . From this prescription, Sb and Sb(0) are given by 

S B = — \ G(q) , S B (0) = — -J—±-G(q) . (2.51) 



3 Results and discussion 

3.1 Aa interaction 

In order to gain the overview of the zero-momentum Wigner transform given 
in Eq. (2.37), we will show in Fig. 1 the central (V#(r, 0)) and LS (V^ s (r, 0)) 
components for the B$a interaction with B$ = N, A, £ and S, when the 
quark-model G-matrix B 8 B 8 interaction by fss2 are employed with the con- 
tinuous choice for intermediate spectra. In this case, we choose qi = and the 




01234 01234 
r (fm) r (fm) 

Fig. 1. The central component of Fig. 2. The same as Fig. 1, but for the 
the zero-momentum Wigner transform LS component Vw LS (r, 0). 
Vw C (r,0) for the B^a Born kernel, cal- 
culated from the quark-model G-matrix 
B$Bs interactions by fss2. The Fermi mo- 
mentum used in the G-matrix calculation 
is k F = 1-35 fm -1 and the continuous 
choice is used for intermediate spectra. 
The (0s) 4 shell-model wave function with 
the h.o. size parameter v = 0.257 fm -2 is 
used for the a cluster. 
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relative momentum q is assigned to the integral mesh point q in Eq. (2.37). The 
Fermi-momentum used for the G- matrix calculation is kp = 1-35 fm -1 , and the 
(Os) 4 shell-model wave function with the h.o. size parameter v = 0.257 fm -2 
is used for the a-cluster folding. We note that the correct treatment of the 
total cm. motion for the B 8 a system is very important, since the cm. cor- 
rection in the standard approach is of the order of 1/A ~ 1/4. As the result, 
the zero-momentum Wigner transform becomes very much short-ranged and 
deep with the interaction range about R = 1.2 A 1 / 3 ~ 2 fm. The a-particle 
density in Eq. (2.44) is more compact than the density including the cm. mo- 
tion, p(r) = 4(2v/tc) 3/2 e~ 2L,r2 , and the central density is about (4/3) 3/2 ~ 1.5 
times larger. The extremely large Na Wigner transform in Fig. 1 is because 
of the factor 2 in Eq. (A. 9), and also because the other nucleon-exchange term 
than the knock-on term and the effect of the exchange normalization kernel 
are neglected. We will deal with this case in a separate paper. 

We list in Table 1 the values of the zero-momentum Wigner transform at the 
origin r — 0, and the Scheerbaum-like factor S B in Eq. (2.47) for all the possi- 
ble combinations of G-matrix calculations for the models, fss2 and FSS, and for 
the QTQ and continuous choices for intermediate spectra. The corresponding 
values for the single-particle potentials, U B (qi), and the Scheerbaum factors, 
Sb(qi), at q± — in symmetric nuclear matter with kp = 1.35 fm -1 are listed 
in Table 2. By comparing these results, we obtain the following findings: 

Table 1 

Values of the zero-momentum Wigner transform at the origin r = 0, V^-(0,0) and 
V^ s (0,0), and the Scheerbaum-like factor, S B , in Eq. (2.47). The unit is in MeV 
for V#(0,0) and V^ 5 (0,0), and in MeVfm 5 for S B . The models are fss2 and FSS, 
and qtq and cont. imply the QTQ and continuous choices for intermediate spectra, 
respectively, used in the G-matrix calculations. 





B 


FSS 




fss2 








qtq 


cont 


qtq 


cont 




N 


-112.8 


-121.1 


-112.6 


-120.0 


t#(0,0) 


A 


-26.01 


-26.03 


-28.09 


-29.31 




S 


59.33 


48.72 


51.76 


47.24 






7.97 


10.73 


20.83 


20.38 




N 


-65.60 


-67.65 


-63.89 


-65.31 




A 


-14.56 


-14.70 


-17.89 


-18.91 




S 


-28.07 


-23.45 


-24.32 


-22.92 






18.59 


25.55 


7.71 


8.59 




N 


-51.00 


-52.24 


-52.15 


-53.49 


S B 


A 


-5.32 


-5.04 


-13.22 


-13.96 




E 


-34.87 


-30.36 


-31.69 


-30.89 






20.15 


26.30 


6.77 


7.62 
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Table 2 

The depth of the single-particle potentials, Ub(qi), and the Scheerbaum factors, 
Sb(qi), at qi = obtained by the angular-averaged G-matrix calculations of the 
quark-model potentials in symmetric nuclear matter. The Fermi-momentum /cp = 
1.35 fm _1 is used. The models are fss2 and FSS, and qtq and cont. imply the QTQ 
and continuous choices for intermediate spectra, respectively. 





D 
JD 


r 3D 




fcc9 
ISSZ 








qtq 


cont 


qtq 


cont 




N 


-79.8 


-89.3 


-80.6 


-88.9 


17b (0) 


A 


-42.9 


-46.3 


-44.8 


-48.4 




£ 


16.1 


17.3 


9.5 


7.3 






-14.9 


-20.8 


-5.3 


-8.0 




TV 


-40.3 


-41.4 


-41.3 


-42.4 


5 B (0) 


A 


-3.9 


-3.6 


-10.0 


-10.6 




E 


-27.4 


-22.6 


-24.1 


-23.2 






14.6 


21.4 


4.6 


5.8 



1. The Aa central Wigner transforms are fairly shallow in comparison with 
U\(0) in the symmetric nuclear matter calculations. Namely, |V^(0,0)| 
is less than 30 MeV in all the cases, while |C/a(0)| is more than 40 MeV. 

2. The Xa and Sa central Wigner transforms are repulsive, although £/g(0) 
is attractive. In particular, V^(0,0) for £a is strongly repulsive, which 
is due to the quark-model prediction of the repulsive SiV(J = 3/2) 3 Si 
interaction. These characters are the result of strong isospin dependence 
of the EiV and SiV interactions, which is discussed in the next subsection. 

3. The LS component, (0, 0), for the Act interaction is by no means 
extremely small, in comparison with that for the £a interaction. The 
ratio is^only 70 - 80 % for fss2 and 50 - 60 % for FSS. On the other 
hand, Sb factors reflect the characteristics of Sb(0); namely, Sb is about 
20 - 30 % larger than £#(0), which is about equal to the enhancement 
factor ^ ~ 1.25 - 1.35 in Eqs. (2.47) and (2.51). 



We find that the zero-momentum Wigner transform is not good enough to 
define phase-shift equivalent local potentials, especially for the Aa interaction, 
In order to see this clearly, we examined the Wigner transform for several 
effective Aa potentials, which can be easily derived from the Aa Born kernels 
given in Appendix B of Ref . [7] . The effective AiV potentials we examined are 
the Sparenberg-Baye potential (SB potential) [7] given by 



van = 



vCE)- 



u 2- 
2 + — 



u 



(3.1) 



with 
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Fig. 3. Solutions of the transcendental 
equation Eq. (3.3), obtained from the Aa 
Wigner transform for the various effec- 
tive AN potentials. SB is the Sparen- 
berg-Baye potential, Eqs. (3.1) and (3.2), 
NSC, ND, NF are the simulated versions 
of the Nijmegen potentials, and JA, JB 
are those of the Jiilich potentials. The en- 
ergy E = —3.12 MeV is assumed. The 
(0s) 4 shell-model wave function with the 



h.o. size parameter v 
used for the a-cluster. 



0.257 fm -2 is 



[ 1 E) = -128.0 exp(-0.8908 r 2 ) 
[ 3 E) = -56.31 / exp(-0.7517 r 2 



> 



-30 



G-matrix k F =1.35fm" 




fss2 (cont) 
fss2 (qtq) 

FSS (cont) 
FSS (qtq) 



2 3 
R (fm) 



Fig. 4. The same as Fig. 3, but 
for the Wigner transform, V^(r,q) in 
Eq. (2.36), calculated from the quark- 
model G-matrix interactions. The Aa 
Born kernels are calculated from the 
quark- model G-matrix BgBg interactions 
by fss2 and FSS. The Fermi momen- 
tum used in the G-matrix calculation is 
kp = 1.35 fm -1 , and qtq stands for the 
QTQ prescription and cont for the con- 
tinuous choice for intermediate spectra. 

- 1015 exp(-5.383 r 
+ 1072 exp 



-13.74 r 2 ) , 



(3.2) 



and the G-matrix simulated AN forces [18] generated from the various OBEP 
potentials, NS (Nijmegen soft-core model NSC89), ND (hard-core model D) , 
NF (hard-core model F), JA (Jiilich model A), and JB (model B). By choos- 
ing the parameters, u = 0.94687 and / = 0.8923 in Eqs. (3.1) and (3.2), 
we can correctly reproduce the A separation energy in ^He: i?g xp (^He) = 
—3.12 ± 0.02 MeV. The strength of the short-range repulsive term (the third 
component) of the NS - JB potentials are slightly modified from the original 
values, in order to reproduce this value. (See Ref. [7].) The phase-shift equiva- 
lent local potential in the semi-classical WKB-RGM approximation [32,33,30] 
is calculated by solving the transcendental equation 



U cS (R) = G W (R, [E - U cS {R)} 



(3.3) 



for some specific energies E, where G (R, q) is assigned to V^}(r, q) in Eq. (2.36) 
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with R = r = \r\ and q = |q|o We here study only the S wave, by neglecting 
the usual semi-classical centrifugal term h 2 (£ + l/2) 2 /2/iA a . The centrifugal 
potentials are included only in the Schrodinger equation. This is a plausible 
approximation, since the LS term of the Aa interaction is very small. The ob- 
tained Aa effective local potentials with E = —3.12 MeV are plotted in Fig. 3 
for SB - JA potentials. The depth of U e s is tabulated in Table 3, together with 
the q 2 value at R — 0, determined self-consistently. The bound-state energy 
of this effective local potential, E B , is calculated by solving the Schrodinger 
equation 

-V—J^ + U ^ V(R) = E B y(R) ■ (3.4) 
{any 

We find that the bound-state energies of SB - JA potentials are too small, 
compared with the exact value Eb (exact), which is obtained by solving the 
Lippmann-Schwinger equation using the Aa Born kernels. Namely, Eb from 
U e s(R) is from 1.3 MeV to 1.7 MeV too small in magnitude, except for the 
rather moderate difference 0.74 MeV in ND. Figure 3 shows that this difference 
is related with the interaction range of U e g(R); i.e., the range of ND is long 
while the others are short. This poor result of the WKB-RGM approximation 
for the Aa interaction is probably related to the very strong nonlocality (or 
momentum dependence) originating from the P r term in Eq. (3.1). In order to 
see this, we artificially changed the Majorana exchange mixture parameter u in 
Eq. (3.1) and compared Eb obtained by the Wigner transform technique and 
by the exact method using the Aa Born kernel. Table 4 shows this comparison. 
The case u = 2 corresponds to pure Wigner-type Aa interaction, which gives 
a local Aa potential and complete agreement between the two methods. Once 
we decrease the u value and introduce the Majorana component, the Wigner 
transform technique loses the attractive effect of nonlocality very much and 
eventually reaches at a very weak effective local potential with no bound state 
before u = (the strength of the odd force =— (the strength of the even force)). 
Our case is just in the middle of these two extremes, which corresponds to 
the approximate Serber-type interaction with a weak odd force. On the other 
hand, the exact solution is almost independent of the u value, which implies 
that the A-particle is bound to the a-cluster in the almost S wave. 

The situation is almost the same even with the Wigner transform approach of 
the quark-model G-matrix interaction. We show in Fig. 4 solutions of the tran- 
scendental equations obtained from the Wigner transform of Aa Born kernels. 
Here we also assumed E = —3.12 MeV. We find that different prescriptions 
for the G-matrix calculations, the QTQ (qtq) or the continuous (cont) choice 
for intermediate spectra, give essentially same result with a slightly smaller 
attraction for the QTQ. The model FSS gives a little weaker attraction than 
fss2. The bound-state energies listed in Table 3 show that the Schrodinger 



2 For the effective AN forces, the approximate formula Eq. (2.36) gives an exact 
result, since V (k, q) is (k • q)-independent. 
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Table 3 

The depth of the effective local potential U e s(0), obtained by solving the transcen- 
dental equation Eq. (3.3) with E = —3.12 MeV. The q 2 value at R = and a special 
value V^(0, 0) of the Aa Wigner transform at r = and q = are also given. The 
eigenvalue Eb is obtained by solving the S'-wave Schrodinger equation Eq. (3.4) for 
U e s(R). The heading Eb (exact) indicates the exact eigenvalue, calculated from the 
Aa Born kernel. The quark-model G-matrix interactions are fss2 and FSS both in 
the QTQ (qtq) and continuous (cont) prescriptions for intermediate spectra. The 
depth of the LS potential U^(0) is calculated from the q 2 value determined by 
using only the central force. 



model 


v£(o,o) 


q 2 


U eS (0) 


ug{0) 


Eb 


Eb (exact) 




(MeV) 


(fm- 2 ) 


(MeV) 


(MeV) 


(MeV) 




SB 


-58.86 


1.062 


-27.15 




-1.86 


-3.12 


NS 


-71.00 


1.269 


-31.85 




-1.74 


-3.12 


ND 


-32.08 


0.613 


-16.98 




-2.38 


-3.12 


NF 


-61.91 


0.959 


-24.83 




-1.78 


-3.12 


JA 


-74.35 


1.168 


-29.56 




-1.60 


-3.12 


JB 


-83.71 


1.312 


-32.81 




-1.43 


-3.12 


fss2 (cont) 


-29.53 


0.731 


-19.66 


-17.22 


-2.94 


-3.62 


fss2 (qtq) 


-28.04 


0.657 


-17.99 


-16.24 


-2.15 


-2.75 


FSS (cont) 


-26.35 


0.673 


-18.34 


-11.50 


-2.62 


-3.18 


FSS (qtq) 


-25.95 


0.607 


-16.87 


-12.30 


-1.78 


-2.29 



Table 4 

Comparison of the S'-wave bound-state energies calculated from the Wigner trans- 
form technique (Eb) and the Lippmann-Schwinger approach of the Aa Born kernel 
(Eb (exact)), when the Majorana exchange parameter u in Eq. (3.1) is changed from 
u = 2 (pure Wigner) to u = (pure Majorana). 



u 


E B 


Eb (exact) 




(MeV) 


(MeV) 


2 


-3.219 


-3.219 


1 


-1.864 


-3.120 


0.5 


-1.272 


-3.079 





unbound 


-3.041 



equation of U e g(R) gives smaller energies by 0.5 ~ 0.7 MeV than the exact 
method using the Aa Born kernel. In fact, the bound-state energy for ^He, 
obtained by solving the Lippmann-Schwinger equation is —3.62 MeV for fss2 
(cont) and —3.18 MeV for FSS (cont). These values are by no means too large, 
compared with the experimental value £^ xp (iHe) = -3.12 ±0.02 MeV. Table 
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Fig. 5. The central component of 
the zero-momentum Wigner transform, 
V^-(R, 0), and the solution of the tran- 
scendental equation, U c s(R), for the 
quark-model G-matrix interaction fss2. 
The continuous choice for intermediate 
spectra is used for the G-matrix calcu- 
lation. 



Fig. 6. The same as Fig. 5, but for the 
LS component. The local momentum q 
is determined self-consistently by using 
only the central components as in Fig. 5. 



3 also shows that the attraction is reduced when the zero-momentum Wigner 
transform is converted to the effective local potential. This is depicted in Fig. 5 
for a typical example of the fss2 prediction with the continuous prescription. 
The small difference of V^(0, 0) values in Table 3 from those in Table 1 is be- 
cause r = 0.1 fm and q = 0.013 fm _1 are actually used in Table 1, instead of 
r = q = 0, and also because a spline interpolation for q is applied to Vo(k, q) 
in Eq. (2.36), to facilitate the Wigner transform for an arbitrary q. The LS 
component U^(R) depicted in Figs. 6 and 7, and the depth U^(0) in Table 3 
are calculated from the LS Wigner transform in Eq. (2.36) by using local mo- 
mentum q determined self-consistently for each R, with respect to the central 
component. First we find in Fig. 6 that the modification of the zero- momentum 
LS Wigner transform to the effective local potential is comparatively small. 
Secondly, Fig. 7 shows that the model FSS gives a shallower LS potential than 
model fss2, which is a result of the strong cancellation between the ordinary 
LS and the antisymmetric LS (LS^) forces in the model FSS. In Fig. 8, we 
compare the central components of the effective local potentials obtained by 
the quark-model G-matrix interaction and by some of the effective AN forces. 
We find that the G-matrix prediction is rather long-ranged and situated in 
the middle of ND and SB predictions. The shape of our Aa central potential 
is rather similar to one of the phenomenological potentials in Ref. [19] (B of 
Fig. 5), and to the effective potential derived from a five-particle microscopic 
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wave function with hard core correlations in Ref. [20] (the dashed curve of 
Fig. 6). It is convenient to parametrize the obtained Aa potentials in simple 
Gaussian functions. For example, the effective local potential at E — —3.12 
MeV, predicted by fss2 (cont) in Fig. 8, is expressed as 

U(R) = -(19.42 + 19.50 i? 2 ) e" ' 5145 ^ (MeV) , (3.5) 

with the bound-state energy E B = —2.95 MeV (which corresponds to —2.94 
MeV in Table 3). 

We list in Table 5 the depths of the effective local potentials U c g(0) for positive 
energies, E cm . = 10, 30 and 50 MeV, and the S'-wave phase shifts S^(E) 
obtained by solving the Schrodinger equation. The exact phase shift S (E) 
(exact) is also shown, which is obtained from the Born kernel. Here we again 
find the effective local potentials are too shallow, especially for the effective 
AN forces. The phase shift difference is about 5 - 11° for the effective AN 
forces, while 3-5° for the quark-model G-matrix interactions. This implies 
that we need a readjustment of the effective local potentials of the order of 10 
MeV, in order to reproduce the correct magnitude of the phase shifts. 

Finally, we will make a brief comment on the choice of the hp value in the 
present framework. If we choose a smaller kp, the A s.p. potential in the 
symmetric nuclear matter becomes shallower and the Aa central interaction 
becomes more attractive. The first feature reduces the magnitude of the start- 




012345 012345 
R (fm) R (fm) 

Fig. 7. The same as Fig. 4 but for the Fig. 8. Comparison of effective local 
LS components. The local momentum is potentials predicted by the quark-model 
determined only by the central force. G-matrix interactions fss2, FSS (with the 

continuous choice for intermediate spec- 
tra) and by some of the effective AA 
forces, SB and ND potentials. 
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Table 5 

The depths of the effective local potentials, U e s(0), and the S-wave phase shifts, 
5^ {E), obtained by solving the Schrodinger equation. The exact phase shift, Sq(E) 
(exact), is also shown, which is obtained by solving the Lippmann-Schwinger equa- 
tion from the Born kernel. 



model 


E 

J — 'cm. 


U pf r(0) 

w en v / 


S^(E) 


Sq(E) (exact) 




(MeV) 


(MeV) 


(dee) 


(deg) 




10 


-20.74 


62.74 


72.20 


SB 


30 


-11.90 


29.01 


36.28 




50 


-4.26 


15.37 


20.69 




10 


-25.16 


61.67 


72.88 


NS 


30 


-15.58 


28.12 


37.10 




50 


-6.74 


13.46 


20.46 




10 


-10.69 


62.14 


68.39 


ND 


30 


-1.78 


26.89 


32.90 




50 


6.29 


12.83 


17.94 


fss2 


10 


-15.33 


70.58 


74.82 


(cont) 


30 


-9.46 


36.24 


40.54 




50 


-4.04 


22.47 


25.75 


fss2 


10 


-13.18 


64.92 


70.07 


(qtq) 


30 


-6.63 


31.28 


35.50 




50 


-1.09 


18.18 


21.54 


FSS 


10 


-14.60 


69.70 


73.49 


(cont) 
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36.88 


40.50 
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-4.96 


23.01 


26.62 


FSS 


10 


-12.21 


62.64 


67.88 


(qtq) 


30 


-6.10 


30.25 


34.30 




50 


-1.27 


18.14 


21.20 



ing energy in the AN G-matrix, resulting in more attractive AN G-matrix. 
Although this change is compensated by the smaller volume in the phase-space 
integral to calculate the A s.p. potential, such a mechanism does not work in 
the present calculation of the Aa interaction. Thus the Aa central interaction 
after the a-cluster folding becomes more attractive for a smaller hp- The same 
situation is observed in Fig. 4 of Ref. [25] for their E(3AQ potential. For exam- 
ple, if we change kp = 1.35 fm" 1 to 1.20 fm _1 (which corresponds to the 70% 
of the normal density), the depth of the A s.p. potential, t/ A (0), in Table 2 is 
reduced from —48.4 MeV to —37.5 MeV for the model fss2 in the continuous 
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choice. On the other hand, the depth of the effective local potential, U e g(0), 
and Eg (exact) in Table 3 changes from —19.66 MeV to —21.74 MeV and from 
—3.62 MeV to —4.54 MeV, respectively!!] This implies that the self-consistent 
mechanism of the starting-energy dependence, which is not properly taken into 
account in this paper, is in fact very important. This finding is in accord with 
the importance of the Brueckner rearrangement effect discussed in Ref. [14]. 
Since the purpose of the present study is not to examine the change of the 
a-cluster, it would be safe to assume the standard value k F = 1.35 fm -1 , in 
order to examine the qualitative features of the B 8 a interaction. 

3.2 Sa and Sa interactions 

As seen from Fig. 1, the Sa and Sa interactions are repulsive. This does not 
mean that the E and H potentials are always repulsive. The structure of spin- 
isospin factors for the Sa and Ea systems is especially simple (see Eq. (2.10)), 

3 For the LS component, the U£g(Q) values in Table 3 change from —17.22 MeV 
to -17.23 MeV for fss2 and from -11.50 MeV to -9.95 MeV for FSS, when k F = 
1.20 fm -1 is used in the continuous choice. 




012345 012345 
R (fm) R (fm) 



Fig. 9. The central components of Fig. 10. The same as Fig. 9, but for fss2. 
the zero-momentum Wigner transform, 
V w C (r,0) with R = \r\, for the Sa Born 
kernel, calculated from the quark-model 
G-matrix B^B^ interactions by FSS. The 
Fermi momentum kp = 1.35 fm -1 and 
the continuous choice are used for the 
G-matrix calculation. The h.o. size pa- 
rameter v = 0.257 fm~ 2 is used for the 
a-cluster. 
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Fig. 11. The same as Fig. 9, but for the Fig. 12. The same as Fig. 11, but for fss2. 
LS components. 

which is due to the spin-isospin saturated character of the a-particle. It is, 
therefore, important to examine each isospin component separately, in order 
to gain some insight to other possibilities of unknown hypernuclei. Here we 
discuss some qualitative features of these interactions, based on the symmetry 
properties of the B 8 B 8 interactions predicted by the quark-model interactions, 
FSS and fss2. When the interaction is repulsive, the transcendental equation 
Eq. (3.3) sometimes does not have its solution, since the square of the local 
momentum, q 2 , becomes negative. Since the extension of the Wigner transform 
Eq. (2.36) to the negative q 2 is not easy numerically, we only discuss the zero- 
momentum Wigner transform in this subsection. 

Figures 9 and 10 show the isospin components with 1 = 1/2 and 1 = 3/2 for 
the Sa central interaction, predicted by FSS and fss2, respectively. We find 
that both interactions have some amount of attraction originating from the 
3 S\ channel of the 1 = 1/2 E7V" interaction. This channel becomes attractive 
due to the very strong AiV-XiV coupling by the one-pion exchange tensor 
force. On the other hand, the 3 Si state of the 1 = 3/2 channel is strongly 
repulsive due to the Pauli principle at the quark level. We find from Figs. 
9 and 10 that this repulsion is so strong that almost no attraction from the 
1=1/2 channel remains in the Sa interaction. 

On the other hand, the two isospin channels, 1=1/2 and 3/2, yield fairly 
large LS force for the Set interaction. This can be seen in Figs. 11 and 12. In 
the isospin 1 = 3/2 channel, 3 Pj states are classified to the flavor symmet- 
ric channel with the SU 3 label (22). It is, therefore, plausible that the same 
mechanism as the NN 1 = 1 channel yields very strong LS force. On the 
other hand, a part of the LS force from the isospin 1 = 1/2 channel is due to 
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Fig. 13. The same as Fig. 9, but for the Fig. 14. The same as Fig. 13, but for fss2. 
Ea interaction. 
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Fig. 15. The same as Fig. 13, but for the Fig. 16. The same as Fig. 15, but for fss2. 
LS components. 

the flavor-exchange process between (11) and (ll) s SU 3 configurations. This 
process is accompanied with the spin flip between S = and 1, and yields 
very strong LS^ force. The LS and LS^ forces reinforce each other in just 
opposite way to the AiV interaction, and the resultant So; LS force becomes 
almost 3/5 of the Na LS force, as seen from Sb in Table 1. This ratio is 
almost the same as that of the Scheerbaum factors Sb = Sb(0) in symmetric 
nuclear matter. (See Table 2.) 



Figures 13 and 14 show the isospin components with 1 = and 1=1 for the 
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Ha central interaction, predicted by FSS and fss2, respectively. Here again we 
find the situation that the attractive nature of the 1 = component is largely 
canceled by the repulsion in the 1 = 1 channel. However, this cancellation is 
not strong especially in the model FSS, and we have almost 5 MeV attraction 
around R = 2 fm. In fss2, the height of the central repulsion in the 1 = 1 
channel is almost 20 MeV, and we can expect a few MeV attraction in the 
surface region. These long-range attractions may have some influence to the 
atomic orbit between H~ and a. The Ha potential rather similar (though more 
attractive) to that by FSS is presented by Myint and Akaishi in Ref. [26]. It 
should be noted that the origin of the repulsion in the 1 = 1 channel is the 
Pauli forbidden state (ll) s in the 1 S'o state and the almost Pauli forbidden 
state (30) in the 3 Si state. However, the coupling with the AS channel is 
very important, which may cause the long-range attraction even in the 1=1 
channel. 

The LS components of the Ha interaction is repulsive, which is clearly seen 
in Figs. 15 ans 16. This LS force is fairly strong, especially for FSS. The 
magnitude is almost 1/3 of the Na LS force with the opposite sign. These 
features are very similar to the LS force in symmetric nuclear matter, as seen 
in Table 2. 



4 Summary 

The SUq quark-model baryon-baryon interaction (fss2, FSS) [1,2,3] is a uni- 
fied model which describes all the baryon-octet baryon-octet (B 8 B$) interac- 
tions in a full coupled- channel formalism. For the nucleon-nucleon (NN) and 
hyperon-nucleon (YN) interactions, all the available experimental data are 
reasonably reproduced. It is, therefore, interesting to study B$a interactions 
in a microscopic framework under a simple assumption of the (0s) 4 harmonic- 
oscillator shell-model wave function for the a-cluster. In this study, we have 
used the result of G-matrix calculations for symmetric nuclear matter, as an 
input for the two-body interactions for the a-cluster folding. Since the resul- 
tant B$a interactions are rather insensitive to the Fermi- momentum &p in 
the G-matrix calculations (except for the AN LS force for FSS), we have as- 
sumed hp = 1-35 fm -1 . The other G-matrix parameters, the center-of-mass 
(cm.) momentum K of two interacting particles and the starting energy u, 
are treated unambiguously in the total cm. frame of the B$a system. This can 
be achieved by considering the transformation of the matrix elements in the 
momentum representation from the initial (qj and final (q^) momenta to the 
momentum transfer (k = qj — q { ) and the local momentum q = (q f + <jfj)/2. 
If one uses this transformation at the level of G-matrix, the procedure of the 
a-cluster folding becomes extremely simple both for the central and LS com- 
ponents. The B$a interaction, V n (k,q) with Q = C, LS, represented by k 
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and q is then transformed back to the Bga Born kernel, V (qj, q { ), through 
the inverse transformation. This procedure is also convenient to calculate the 
Wigner transform, V^(r, q), which is simply a Fourier transform of V n (k, q). 
By solving the transcendental equation for V$(r, q), we can obtain an energy- 
dependent local potential of the B 8 a system in the WKB-RGM approximation 
[32,33,30]. 

In this paper, we have applied the present formalism to the Aa, Ea and Ea 
systems. Applications to the Aa system will be discussed in a separate paper, 
since this system involves an extra nucleon-exchange term, in addition to the 
direct and knock-on terms. In the Aa system, the WKB-RGM approximation 
is rather poor due to the very strong momentum dependence of the exchange 
knock-on term. This term appears even for the effective AA force, which is 
traced back to the strangeness exchange processes. We find that the Aa central 
potentials predicted by various quark-model G-matrices are very similar to 
each other, irrespective of a specific model, fss2 or FSS, and the QTQ or 
continuous choice for intermediate spectra. At the bound-state energy, E = 
—3.12 MeV, they are long-range local potentials with the wine-bottle shape, 
having the depth less than 30 MeV. They are very similar to the Aa potential 
obtained from the effective AA potential ND [18], simulating the Nijmegen 
hard-core model D. The Aa bound-state energies are calculated by solving 
the Lippmann-Schwinger equation of the Aa Born kernel. These are —3.62 
MeV for fss2 (cont) and -3.18 MeV for FSS (cont), when k F = 1.35 far 1 
is used. It should be noted that the depth of the single-particle potential for 
A in symmetric nuclear matter is 48 MeV for fss2 (cont) and 46 MeV for 
FSS (cont) [1]. The fact that the predicted Aa bound-state energies are by 
no means too large, in comparison with the experimental value E^ p (jjie) = 
—3.12 ± 0.02 MeV, implies that the proper treatment of the cm. motion of 
the Aa system is very important. It is also important to note that, in the 
present G-matrix approach, the AN-Y>N coupling by the very strong one- 
pion exchange tensor force is explicitly treated, the lack of which is known to 
lead to the so-called overbinding problem of the Aa bound state. The energy 
loss predicted by the a-cluster rearrangement effect [14] through the starting- 
energy dependence of the G-matrix needs further detailed analyses. On the 
other hand, the Aa LS potentials are rather model dependent. In the model 
fss2, the depth of the LS potential is —16 ~ —17 MeV, while in FSS about 
—12 MeV. The interaction range of the FSS LS potential is also very short, 
which leads to a small Scheerbaum-like factor S\ for the Aa LS force. We will 
show in a separate paper, the strength of the LS force is further reduced for 
smaller values of k F , if FSS is used. The very small spin-orbit splitting of the 
^Be [34,35] can be reproduced in the aaA Faddeev calculation, using the aa 
RGM kernel and the Aa LS Born kernel predicted from the FSS G-matrix 
with kp = 1.25 fm _1 . 

Based on the reasonable reproduction of the Aa interaction properties, we have 
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examined the real parts of the Ea and Ha interactions in the Wigner transform 
technique. Since these interaction are repulsive, we have examined only the 
zero-momentum Wigner transform as the first step. In the Ea interaction, the 
attractive effect from the isospin 1 = 1/2 EA channel is completely cancelled 
out by the repulsion from the 1 = 3/2 EA channel. The origin of this strong 
repulsion is the UN (I = 3/2) 3 Si channel, which contains the almost Pauli- 
forbidden SU3 (30) component for the most compact (Os) 6 configuration. On 
the other hand, the Ha interaction is less repulsive because of the appreciable 
attraction originating from the 1 = EN channels. The Ha zero-momentum 
Wigner transform predicted by FSS yields about —5 MeV attraction around 
R = 2 fm, while the attraction of fss2 is about —3 MeV. These long-range 
attractions may have some relevance to the formation of atomic bound states 
for the H~a system. As to the spin-orbit interaction, the two isospin channels 
of the EA interaction give fairly strong attractive Ea LS forces, yielding 
almost 3/5 of the Na LS force. On the other hand, Ha LS force is repulsive 
and the magnitude is 1/8 ~ 1/2 of the Na LS force. We will show in the 
next paper, the present Na LS force is consistent with the observed P-wave 
splitting of the 3/2" and l/2~ excited states of 5 He. 

It should be noted that the overall repulsive character of the Ea and Ha in- 
teractions is related to the spin-isospin saturated character of the a-cluster. 
The strong isospin dependence of the EA and HA interactions, namely, re- 
pulsive for the EJV(7 = 3/2) and EN (I = 1) channels and attractive for the 
EA(7 = 1/2) and EN (I = 0) channels, leads to a possibility of attractive 
features in some particular spin-isospin channels for systems of E, H and the 
s-shell clusters [23,24,25]. One of the examples of such systems is the isospin 
1 = 1/2 and spin 5 = state of f^rle, in which the strong repulsion of the 
UN (I = 3/2) 3 Si channel does not contribute [36] and a quasi-bound state is 
in fact observed experimentally [37,38]. Applications of the present approach 
to such systems are under way. 
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A Invariant G-matrix for the most general B 8 B 8 interaction 

In this Appendix, we will define the invariant G-matrix for the most general 
B 8 B 8 interaction. The B 8 B 8 channels in the bra side (7) and ket side (a) are 
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specified by [30] 



7 : 



i(ll) Cl i(ll)c 2 



SS Z YII Z ;V , 

S'S g 'YII g -V , (A.l) 
in the isospin basis. For example, the spin-flavor functions are given by 



a = 



^(11)01^(11)02 



= Xss z [BiB 2 ]u z , 
[B 1 B 2 ]f Iz 



y/2(l + 6 CltC2 ) 



{[B 1 B 2 \ n .+V{-\)™*- I [B 2 B 1 \ IIm } , (A.2) 



where the isospin-coupled flavor wave functions, [B 1 B 2 \u z , are lexicograph- 
ically ordered and the first baryon is numbered always 1 and the second 2 
(i.e., [B 2 B 1 ] IIz = [S 2 (l)Si(2)] 77 J. The flavor symmetry basis, [B l B 2 ] 1 ? Iz with 
V — 1 and —1, gives 



\ B i B 2\ IIz 
[ B 3 B i\n z 



'1 + 5, 



C1,C2 



'1 + Sa 



E 

7" 



4J// Z 



(A.3) 



The matrix element of the G-matrix in the isospin basis and its partial-wave 
decomposition are given by [31] 



G ia (p, pf; K, u) = ( X ss z [B.B^jGip, p'; K, cu)\ X s's z > [ B z B $'i z ) 

= E ^G J ^ aS , Sz/ (p,p'-,K,uj)J2^rnSS z \JM)Y em (p) 

JMW m 

xJ2(Z'rn'S'S z '\JM)Y* m ,(p') . (A.4) 

m' 

Here the prime symbol on implies that the summation is only for such 
quantum numbers that satisfy the generalized Pauli principle: 



iy- s v = (-i) t (-iy- li 'v' = -i . 



1-S't>/ 



(A.5) 



We multiply Eq. (A.4) with \xss z ) and (xs'S z '\ from the left- and right-hand 
sides, respectively, and take a sum over SS Z and S'S Z . Then we find 
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([B 1 B 2 ^ Iz \G(p,p';K,io)\[B 3 B 4 }Z z ) 

= E ^ G is£, a s's z '(P^P'i K ^)y(es)JM(p]spm)yi e , sf)JM (p';spm) , 

JMU'SS 1 

(A.6) 

where ^5)jm(p; spin) = \¥z(p)xs\.jm is the angular-spin function for the two- 
baryon system. 

Let us denote the two-baryon channels with c = (01,02), a = (01,02), and 
consider 

G ca (p, p';K,uj) 

= ([B 1 B 2 ] nz \G(p,p , -,K,u;)-G(p,-p'-,K,u;)P a P F \[B 3 B 4 } Ih ) . (A.7) 
We can rewrite this by using Eqs. (A. 3) and (A.6): 

G ca (p, p';K,u) 

= y^ + S cl , c ,)(l + S ai , a2 ) E {{[B l B 2 \l}G{p,p'-K,u J )\[B,B^' Iz ) 

-([B^JG^-p'^K^^PFllBzB^i)} . (A.8) 

Here the second term in the brackets gives the same contribution as the first 
term, due to the condition Eq. (A. 5). Thus we find that Eq. (A.7) is actually 
invariant G-matrix, which allows the expressions 

G ca (p,p'; K,u) 

= {[ B^j^G^p'-K^ ) - G(p,-p';K,Lu)P a P F \[B 3 B 4 ] IIz ) 

= V / ( 1 + ^ C2 )( 1 +<W 2 ) E ([B^JG&p'; K,u)\[BzB$' Iz ) 

v,v 

I 

= V / (l + 5 CllC2 )(l+5a 1 ,a 2 ) E 4ff G isi,aS>S z > V' \ K, u) 

JMWSS'VV 

xy(es)jM(P; spin) y^s')jM(P'^ s P in ) 
= \f(l + o" CljC2 )(l +5oi,o 2 ) [go + 9ss (cri • <r 2 ) + h in • (cr x + cr 2 ) 

+h-in- (<ri -<r 2 ) H ] . (A.9) 

The eight independent invariant functions, go, g ss , etc., are expressed by some 
combinations of the partial-wave components of the G-matrix, 
G ise aS'S z '(PiP'i Ki^)-* which are explicitly given in Appendix D of Ref. [31]. 
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